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ABSTRACT 

We present a numerical investigation of the contribution of the presupernova ejecta of Wolf-Rayet stars to 
the environment surrounding gamma-ray bursts (GRBs), and describe how this external matter can affect the 
observable afterglow characteristics. An implicit hydrodynamic calculation for massive stellar evolution is used 
here to provide the inner boundary conditions for an explicit hydrodynamical code to model the circumstellar gas 
psj ■ dynamics. The resulting properties of the circumstellar medium are then used to calculate the deceleration of a 

>^ relativistic, gas-dynamic jet and the corresponding afterglow light curve produced as the shock wave propagates 

through the shocked-wind medium. We find that variations in the stellar wind drive instabilities that may produce 
radial filaments in the shocked-wind region. These comet-like tails of clumps could give rise to strong temporal 
^5 | variations in the early afterglow lightcurve. Afterglows may be expected to differ widely among themselves, 

, depending on the angular anisotropy of the jet and the properties of the stellar progenitor; a wide diversity of 

' behaviors may be the rule, rather than the exception. 

Subject headings: gamma-rays: bursts — ISM: jets and outflows — radiation mechanisms: non thermal — 
' polarization — relativity — shock waves 

[ 1. INTRODUCTION 

■^j- , Over the past six years evidence has mounted that long-duration (> 2 s) gamma-ray bursts (GRBs) signal the collapse of massive 
(N) ■ stars in our Universe. This evidence was originally based on the probable association of the unusual GRB 980425 with a type Ib/c 
i-H supernova (SN; Galama et al. 1998) but now includes the association of GRBs with regions of massive star formation in distant 
galaxies (Paczyhski 1998; Wijers et al. 1998; Fruchter et al. 1999; Djorgovski et al. 2001; Trentham et al. 2002), the appearance 
' of supernova-like bumps in the optical afterglow light curves of several bursts (Bloom et al. 1999; Zeh, Klose & Hartmann 2004 
^ (— I and references therein), lines of freshly synthesized elements in the spectra of a few X-ray afterglows (Piro et al. 2000; Ballantyne 
O ,■ & Ramirez-Ruiz 2001; Reeves et al. 2002), and the first convincing spectroscopic evidence that a very energetic supernova (a 
I ' hypernova) was temporally and spatially coincident with a GRB (Hjorth et al. 2003; Stanek et al. 2003). These observations support 
Q • the idea that long-duration GRBs are associated with the deaths of massive stars, presumably arising from core collapse (Woosley 
■ 1993, Zhang et al. 2003). 

An implication of a massive star progenitor is that the circumburst environment is determined by the mass-loss wind from the 
• • . star. Much of our effort in this paper will therefore be dedicated to determining the state of the circumburst material in various 
. £h ' tyP es °f progenitor scenarios, and describing how this external matter can affect GRB jets propagating through it. It should be noted 
k>( \ here that the afterglows sample a region ~ 10 17 cm in size. Because massive stars are expected to have their close-in surroundings 
'j_J ■ modified by the progenitor winds, we consider both free winds and shocked winds as possible surrounding media for the afterglow 
C$ stage. Detailed hydrodynamic simulations of this interaction are presented in §4. Our computations allow us, for the first time, to 
study the behaviour of circumstellar gas very close to the progenitor star. An understanding of the evolution of a gas-dynamic jet can 
come only through a knowledge of the properties of the medium which it propagates. Calculations of the evolution of a relativistic, 
gas-dynamic jet and expected emission properties are discussed in §5. For completeness, the interactions with either a free wind or a 
constant density medium, as well as the termination shock wave which marks the transition between these two media, are discussed 
in §2 and §3. The role of binarity is briefly addressed in §4.1. We then in §6 discuss the possible variety of afterglow variability that 
is expected from GRB jets expanding in a medium that may be inhomogeneous because of clumping in the Wolf-Rayet (WR) star 
wind. Discussion and conclusions are presented in §7. 



2. THE CIRCUMBURST MEDIUM 

If the progenitors are massive stars then there is an analogy to the explosions of core collapse supernovae, for which there is 
abundant evidence that they interact with the winds from the progenitor stars. In most supernova cases, the radial range that is 
observed is only out to a few pc, such that the mass loss characteristics have not changed significantly during the time that mass 
is supplied to the wind (Chevalier & Li 2000). The density in the wind depends on the type of progenitor. Red supergiant stars, 
which are thought to be the progenitors of Type II supernovae, have slow dense winds. Wolf-Rayet stars, which are believed to be 
the progenitors of Type Ib/c supernovae and possibly of GRBs (e.g. MacFadyen & Woosley 1999), have faster, lower-density winds. 
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The winds from typical red supergiants are slow-moving and dense, with velocities v w w 10-20 km s" 1 and mass loss rates between 
10" 6 and 10" 4 M yr _1 (e.g. Fransson et al. 1996). The winds from WRs, on the other hand, are characterized by mass-loss rates 
M w 10" 5 - 10" 4 M Q yr _1 and velocities v w w 1000-2500 kms" 1 (e.g. Chiosi & Maeder 1986). In a steady, spherically symmetric 
wind, the electron density is 

Mr) = 4^v X, m M 3 X lO^m-V^v^M; 1 , (1) 

where /j, e is the molecular weight per electron and [i e ~ 2 in a helium gas. Here v w = 10 3 v Wi 3kms -1 , r = 10 15 ri5 cm, and M = 
10- 4 M_ 4 M Q yr 1 . 

For this discussion we shall assume the blast wave is adiabatic, i.e. its energy is constant with time, and effectively spherical. This 
means the E here is the isotropic equivalent energy as, for example, derived from the gamma-ray output. Deceleration due to the 
stellar wind starts in earnest when about half the initial energy is transferred to the shocked matter, i.e. when it has swept up 7" 1 
times its own rest mass. The typical mass where this happens is 

M dec = « 5 x lO- 6 £ 5 372" 2 ^0- (2) 

The relativistic expansion is then gradually slowed down, and the blast wave evolves in a self-similar manner with a power-law 
lightcurve. This phase ends when so much mass shares the energy that the Lorentz factor, 7, drops to 1. Obviously, this happens 
when a mass E/c 2 has been swept up. This sets a non-relativistic mass scale 

M NR =^«6x 1O- 2 £ 53 M . (3) 

Beyond this point, the event slowly changes into a classical Sedov-Taylor supernova remnant evolution, leading to a steeper decline 
in the lightcurve (Waxman et al. 1998). 

In the unshocked wind, the mass within radius r is Mr/v w , which combined with equation (2) gives the blast wave deceleration 
radius in a stellar wind: 

r d = ~ 2 x 10 16 E 53 v w . 3 AQ 1 ?cm, (4) 

where £53 = (£/10 53 ) ergs. By contrast, the well-known expression for the uniform-medium deceleration radius is 

r d « 10 18 (£ 53 /» ism ,o) 1/3 72 8/3 cm, (5) 

and so a blast wave in a wind decelerates at a much smaller radius. The deceleration time is given by = r<i/(2j 2 c) in both cases, 
and thus is correspondingly smaller for the wind case. However, because the Lorentz factor decreases as M~ x l 2 (with M being the 
swept-up mass) beyond this point, the mass-starved wind blast wave decelerates much more slowly, and therefore begins to catch up. 
There is therefore not much difference in size between wind and uniform CM blast waves in the most commonly observed interval 
from 0.3 to 10 days after the burst (Fig. 1). Usually in the study of afterglows, one considers either the uniform ambient medium 
case or the 1 /r 2 wind case on its own 6 . However, since the wind of a star meets the interstellar medium (ISM) at some point, the 
density structure is more complex, and it is to this problem that we now turn our attention. 

3. WIND-ISM INTERACTION 

During the evolution of a wind-driven circumstellar shell the system has a four-zone structure (analogous to that of a supernova 
shell; Woltjer 1972). From the inside to the outside these zones are: (a) a supersonic stellar wind with density p(r) = M/Aitr 2 v w \ (b) a 
hot, almost isobaric region consisting of shocked stellar wind mixed with a small fraction of the swept-up interstellar gas; (c) a thin, 
dense, cold shell containing most of the swept-up interstellar gas; (d) ambient interstellar gas of number density n; sm (Fig. 2). 

The wind initially expands unopposed into the ISM with a velocity of about v w <~ Voo, the escape velocity at the sonic point. The 
free expansion phase is considered to be terminated at a time t w ^, when the swept-up mass of the interstellar medium is comparable 
to the mass in the wind. The mass lost by the star is Mt w ^ and the swept-up mass is ^j-(v w t Wt d) 3 ni sm m p fi e . These two masses are equal 
when f Wj d = [3M/(47rv 3 t ,ni sm m / ,/i e )] 1//2 , which is about 100 years for a typical WR wind expanding into an homogeneous ISM. The 
free-expansion phase takes place at the early stages of the evolution of the hot star and occupies a minimal fraction of its lifetime. 
During this time both M and v w are approximately constant, and the wind bubble has reached a radius of ~ 9 x 10 11 M- 5 nJ 2 j vj( 2 cm, 

where M is the mass loss rate in units of solar masses per year, and n; sm is the density of the surrounding medium in units of cm" 3 . 

When the free-expansion phase (a) has ended, the wind encounters an inward facing shock. Kinetic energy is deposited in the 
shocked wind region in the form of heat, 

±^(Av w ) 2 =1.4xl0 5 (^ 
16 k 100 kms" 

6 Although the interstellar and wind models are the two main types of environments considered for afterglows, there is a different scenario involving a massive star in 
which the supernova explosion occurs before the GRB (Vietri & Stella 1998). The supernova would expand into the progenitor wind, creating a complex circumburst 
region in the inner part of the wind. Konigl & Granot (2002) have recently shown, for the case of a pulsar-wind bubble, that the shocked wind has a roughly uniform 
density, similar to that found in the normal interstellar medium. 



Tshock = — ^P(Av w ) 2 = 1.4 x 10 5 ( ,_, - _, ) 2 K, (6) 
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where Av w is the relative speed of the material approaching the shock. Thus, a jump in velocity by 800 km s -1 , which is still well 
below typical terminal wind speeds, produces a 10 7 K gas in the shocked wind region. During phase (b), the material is so hot that 
it causes the contact surface to expand outward more slowly than it would in a freely expanding wind. The ISM that enters the 
outward facing shock is heated to a temperature below 10 6 K, emission of line radiation becomes the dominant cooling process and 
the swept-up gas cools quickly to temperatures of about 10 4 K that can be maintained by the radiation field of the star. The duration 
of the adiabatic expansion phase can thus be estimated by finding the time it takes the expanding gas to cool from r s h OC k ~ 10 7 K to 
10 6 K. Using equation (6), we find that a change in temperature from 10 7 K to 10 6 K corresponds to a change in jump velocity by a 
factor of \/l0. This change in jump velocity corresponds to a phase (b):(a) age ratio of about 6 (Castor et al. 1975). Thus, the age of 
the adiabatic phase is less than about 1000 years. 

The mass of the swept-up material is much larger than that in the hot wind and, because it is cool, it lies in a compressed region. 
Phase (c) persists for as long as the star is able to sustain a powerful wind. The dominant energy loss of region (b) is work against 
the compressed region (c). The compressed region (c) expands because its gas pressure is higher than that of the surrounding ISM. 
Therefore, the expansion is described by the momentum equation, 

^ r [M s (t)v(t)]=4irr 2 (t)P i , (7) 
at 

where P\ is the internal pressure of the compressed region, assuming that most of the swept-up interstellar mass remains in the thin 
shell. M$(t) is the mass of the shell of swept-up material, given by Ms(t) = (4/3)7rr 3 (?)pi. P; is determined by the gas pressure of the 
high temperature gas in region (b). The wind material that enters the backward-facing shock is hot, but the material that enters the 
forward-facing shock is cool. The cooled swept-up material is driven outward by the high gas pressure of the hot bubble. The stellar 
wind adds energy to region (b) at a rate 

L(t) = ^M(t)vlit). (8) 

The internal energy in the bubble is given by the product of the energy per unit mass of the material, (3/2)nkT/ p\ = (3/2)Pi/p\, and 
the total mass of the bubble, (4/3)wr 3 pi. Since the total internal energy of the bubble comes from the energy of the wind, we find 
/* = /.„, \2-rU,)\. 

The expansion of the bubble during the adiabatic phase can be found numerically by using this expression in the momentum 
equation. If the wind power L is roughly constant for a period of time, t, one can write f; = Lt/(2irr 3 ). The resulting solution of 
equation (7) gives r(t) oc f 3//5 . This shows that the shell expands more slowly than would a freely expanding wind. Since the gas 
in the cavity moves subsonically its pressure keeps it approximately at uniform density. The bubble could continue to expand until 
stalled by the pressure of the ISM (Garcfa-Segura & Franco 1996). 

The radius of the wind termination shock at the inner edge of the wind bubble can be found by balancing the wind ram pressure 
with the post-shock cavity pressure. For a star that loses mass at a rate lO^M-eA^yr" 1 with a wind velocity 10 3 v w ,3 km s" 1 in 
interstellar gas with density 10 3 nj sm .3cm~ 3 , we have a inner termination shock radius 

rm = 0AM%\^n£X /5 PC, (9) 
where 10 6 fe is the lifetime of the star in Myr. The density in the uniform shocked wind region, n sw , at late times is given by 

^ M nn£ n>2/5 3/5 -6/5-4/5 -3 /1m 

" sw ~ ^,2„ ~ = 0.06M_£ n-J m3 v w j t 6 cm , (10) 

which shows that even if the progenitor star is embedded in a dense molecular cloud the observed blast wave can propagate in a 
low-density, uniform medium (Wijers 2001). The mass within the 1 jr 1 wind, M t , is 

.. Mr, ,n-4»>13/10 -9/10 -3/10 2/5,. 

M t = — = 3 x 10 4 M 6 vj nj^ t 6 < Mq. (11) 

v w 

Comparison with estimates given in Soderberg & Ramirez-Ruiz (2002) show that if the wind is particularly weak (i.e. M < 
10~ 6 MQyr _1 ) or the surrounding density is high («i sm > 10 3 cm" 3 ), r t falls within the range of the relativistic expansion. Models 
and observations of Galactic Wolf-Rayet stars, however, show that the swept-up shell of a red supergiant material at the outer radius 
is at a distance > 3 pc from the star (Garcfa-Segura et al. 1996a). This radius is sufficiently large that the interaction with the free 
1 jr 1 wind is expected over the typical period of observation of afterglows. 

Among the afterglows that can be interpreted as interaction with a undisturbed stellar wind, the highest density objects are com- 
patible with expectations for the wind from a typical Wolf-Rayet star (Panaitescu & Kumar 2002), but the lower densities imply a 
wind densities that are lower by a factor of 10- 10 2 . As proposed by Wijers (2001), the low-density requirement may be alleviated 
by appealing to a shocked wind, but the observable afterglow transitions due to the blast wave traversing the wind termination shock 
wave (Fig. 2) have not been seen in any afterglow. 

Depending upon the wind history of a massive star during its last few centuries, however, the density structure in this inner region 
(i.e. relevant to the afterglow phase) could be quite complicated as the star enters advanced burning stages unlike those in any Wolf- 
Rayet star observed so far. The non-steady nature of the winds in massive stars therefore leaves open the possibility of interaction 
with denser material at early times (Chevalier et al. 2004). This encourages us to present a detailed account of the large and small 
scale features that may be present in the circumstellar gas distribution. These features, as we argue, can result naturally when one 
takes into account the complete mass-loss history of a massive star. 
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4. WIND-WIND INTERACTION 

The detailed dynamical evolution of the circumstellar material (CSM) around massive stars is complex. Some stages of it do not 
involve major hydrodynamical instabilities and can thus be studied analytically by means of self-similar solutions (Weaver et al. 
1977; Ostriker & McKee 1988; Chevalier & Liang 1989 ; Garcfa-Segura & Mac Low 1995a). However, the frequent occurrence 
of instabilities requires two - or three - dimensional hydrodynamic calculations in order to follow the non-linear evolution of the 
resulting structures (Franco et al. 1991; Blondin & Lundquist 1993; Garcfa-Segura & Mac Low 1995b). 

In an effort to achieve full consistency between stellar and circumstellar evolution, we have performed several computations where 
the time-dependent input for the calculation of the circumstellar gas dynamics is derived from the output of a stellar evolution code 
developed at the University of Gottingen (Langer et al. 1988, Langer 1991). Here we carry out computational simulations with the 
hydrocode ZEUS-3D (version 3.4) developed by M. L. Norman and the Laboratory for Computational Astrophysics. ZEUS-3D is 
a finite-difference, fully explicit, Eulerian code (Clarke 1996) descended from the code described by Stone & Norman (1992). We 
used spherical coordinates for our simulations, with a symmetry axis at the pole, and reflecting boundary conditions at the equator 
and the polar axis. The reader is referred to Garcfa-Segura et al. (1996a,b) for a review of the applied computational methods and 
techniques. 

The preburst stellar wind depends on the evolutionary stages prior to (and during) the Wolf-Rayet stage. For Galactic stars, a 
standard evolutionary track is to start as an O star, evolve through a red supergiant (RSG) phase or luminous blue variable (LBV) 
phase with considerable mass loss, and ending as a Wolf-Rayet star (Garcfa-Segura et al. 1996a,b). At low metallicity, the RSG 
phase may be absent; this may also be the case for some binary stars. As a first example, we follow the dynamics of the circumstellar 
medium around a 29 M Q star at the ZAMS, which evolves (at solar metallicity) through a long-lived RSG stage with prominent 
consequences for the evolution of the circumstellar matter. The 29 M Q stellar model has steady winds during the main-sequence 
(MS) and RSG stages. For that reason, the CSM evolution during these stages is computed in one dimension for a homogeneous 
ISM. For the calculation of the post-RSG evolution the variables such as the temperatures and density are extrapolated onto a two 
dimensional grid. Note, however, that the calculations start before the RSG phase has ended in order to be sure that the RSG shell is 
dynamically stable. The total amount of mass lost by the RSG in the wind is M rsg = 10M Q , with a mass loss rate M rsg = 6 x 10~ 5 MQyr _1 
and wind velocity v rsg ^15 km s" 1 . 

When the fast WR wind v w ~ 3500 km s" 1 starts blowing, it sweeps up the RSG wind material into a shell, which we will refer to 
as the WR shell. The properties of this shell turn out to be very sensitive to the characteristics of the RSG wind. Since the density 
of the RSG wind depends on its velocity as p(r) = M r sg/(47rr 2 v rsg ), a low RSG wind velocity implies a much higher density. The 
expansion of the WR shell is faster for lower RSG wind densities and also for higher wind velocities. The termination shock of the 
WR wind is located at a much smaller radius than that of the hot MS bubble described in the previous section, so we simply assume 
free outflow as the outer boundary condition in order to calculate its dynamical behaviour (Garcfa-Segura et al. 1996a). 

During the early stages, the swept-up WR shell is dense enough to be fully radiative. This ceases to be true, however, after the 
shell has extended to a radius of more than 0.1 pc, since its density decreases so much that its evolution becomes almost adiabatic. 
Correspondingly, the shell is initially thin and therefore subject to Vishniac instabilities (Vishniac 1983); but the increase of shell 
thickness soon suppresses the growth of the instability. This instability operates in dense winds where the cooling is effective enough 
to produce radiative terminal shocks. The RSG wind shell is accelerated by the (less dense) shocked WR wind and is therefore also 
subject to Rayleigh-Taylor (RT) and Vishniac instabilities. In this case we have a radiatively-cooled, thin, accelerated shell and the 
two instabilities must be coupled. The resulting structure after 8,000 yr of evolution since the onset of the WR phase is summarized 
in Fig. 3. 

It is also possible for a fast WR wind to collide with early ejecta before it is decelerated by the RSG shell (or a LBV shell depending 
on the evolutionary pathway of the massive progenitor). To illustrate this, we study the last (unsteady) WR stage of a 60 M Q star. 
Mass loss lowers the stellar mass from 60 M Q all the way down to <~ 4M Q at the time of the final supernova (or GRB). The total mass 
loss includes ~ 32 M lost in the MS and pre-LBV stages, <~ 8 M Q lost in the LBV stage, and <~ 16 M Q in the WR stage. The MS 
and LBV-WR stages have been modeled in a self-consistent fashion by Garcfa-Segura et al. (1996b). These stages have important 
implications for the circumstellar gas at r 3> 0.4 pc. Here, we study the interaction between the fast-moving WR wind (v wr <~ 4000 
km s -1 ) just before the final supernova and the early, slower moving (v wr ~ 1500 km s" 1 ), ejecta 7 . 

The resultant shell will be pushed outward by the central star wind and retarded by the early ejected wind, quickly reaching a 
constant velocity, but increasing in mass. The resulting expansion law for the shell follows from the momentum balance between 
the two components. The analytical theory of Kwok et al. (1978) can give some insight into the problem. If M tj is the mass of the 
resulting shell when it is at a radial distance r, 7 (f), then 

r r 'j(i) m ■ f n+Vit M- 
M u (t)= / — -dr+ / — dr, (12) 

where Vj < v,. If v, ; (f) is the velocity of the shell, then, assuming a completely inelastic collision, the equation of motion may be 
written 

M,fl)*£ = ^i [Vi - Vij f-Mi [Vj - Vij f. (13) 
at Vi v i 

Numerical integration of equation (13), with a substitution for Mij(t) from equation (12) gives the resulting expansion law for the 
shell. The thickness of the shell Ar t j may be found by requiring its internal pressure to balance the pressure from the wind. 

7 The reader is refer to Fig. 2 of Garcfa-Segura et al. (1996b) for the stellar mass-loss rates and wind velocities as a function of time for the 60 Mq model. 
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Such large expansion velocities Av w produce high temperatures (see equation 6) in the post-shock region behind the outer shock 
of the merged shell. This high post-shock temperature and the low density (i.e. low M wr and high v wr ) of the WR wind result in an 
almost adiabatic, hot, low-density WR shell. 

Figure 4 shows the morphology of the smooth WR ring under the assumption that the central star experiences non-spherical mass 
loss close to critical rotation: in other words, a scenario in which a slower and denser wind is confined to the equatorial plane. 
To compute the latitudinal dependence of the wind properties of a star close to critical rotation ideally requires multi-dimensional 
models of the star and its outflowing atmosphere, which are not available. Langer (1998), however, argued that the stellar flux and 
the radius might still vary only weakly from pole to the equator in very luminous stars. We therefore applied equations similar to 
those found by Bjorkman & Cassinelli (1993) for winds of rotating stars in the limit of large distance from the star: 

Voo(0) = Cv«c(l-ft8iii0)* , (14) 
where we set the parameters defined in Bjorkman & Cassinelli (1993) to £ = 1 ,<p = 0.35, = v rot /v cr j t , and v cr ; t = v esc /\/2 = 
[GM^I-k)*,] 1 ' 2 , with M* and being mass and radius of the star, and k standing for the ratio L/L E dd of stellar to Edding- 
ton luminosity (Langer et al. 1998). Under the above conditions, the wind expands more quickly and easily into the lower density 
wind at the poles, while stellar rotation concentrates it toward the equatorial plane, producing a double-lobed structure. In section §5 
we aim to examine the interaction of the relativistic blast wave with these more realistic density wind profiles. 



4. 1 . The Role ofBinarity 

One of the most important questions relating to WR stars is whether they are all members of binary systems or rather do some 
truly single objects exist. Possibly a large fraction of stars are members of binary or multiple systems but most are sufficiently far 
from their companions that their evolution proceeds essentially as if they were single stars. However there is an important minority 
of stars which are close enough that their evolution is dramatically changed by the presence of a companion. It is extremely unlikely 
that the progenitors of GRBs are just very massive, single WR stars. Special circumstances are almost certainly needed. The most 
promising of these is rotation: a rapidly rotating core is the essential ingredient in the "collapsar" model for GRBs (Woosley 1993; 
MacFadyen & Woosley 1999). Massive stars are generally rapid rotators on the main sequence. However, there are many well- 
established mechanisms by which they can lose their angular momentum during their evolution. Thus, it is not at all clear whether 
the cores of massive single stars will ever be rotating rapidly at the time of explosion (e.g. Heger et al. 2005). The simplest way to 
generate the required angular momentum is for the core to reside in a tight binary and be in corotation with the binary. Izzard et al. 
(2004), using detailed binary population calculations, concluded that there are enough binaries where tidal locking could account for 
the observed GRB rates. Most companions, being low-mass MS stars, have small mass-loss rates and are thus unlikely to enhance 
significantly the density around the GRB progenitor. The density structure seen around WR 147 (Contreras & Rodriguez 1999) could 
be a hint of what exists in some GRB progenitors. Overall, we do not expect binarity to significant alter the circumstellar gas density 
profiles along the rotation axis of the collapsing stellar core. There is, however, a small, interesting subclass of binaries in which the 
companion star may be massive (see e.g. Izzard et al. 2004). In this case, the WR star mass loss may also be influenced by the winds 
of the close binary companion. Detailed computations are needed to show the magnitude of the effect. 

The requirement that the total core angular momentum exceed the maximum angular momentum of a Kerr black hole places inter- 
esting limits on the binary period (Izzard et al. 2004; Podsiadlowski et al. 2004). We make the reasonable assumption that the binary 
is circular, and model the core as an n = 3 polytrope. then the binary period must be smaller than P or b ~ 4(M core /2M0)" 1 (/? core /lO 10 cm) 2 
h. This orbit could be tight enough that the core may in fact have been stripped of its helium in a common envelope to form a CO 
core. Alternatively, it may sometimes happen that the core of a very massive star retains the required angular momentum as its outer 
hydrogen layers are blown off in a stellar wind. Indeed, Brown & Bethe (1994) have noted that the cores of stars more massive than 
~ 20M Q will undergo a delayed collapse to form a black hole if the nuclear equation of state is soft. Prompt formation of a black hole 
introduces a mechanism for failure of a core collapse supernova, if success of the shock depends on delayed neutrino heating. The 
corresponding rate of Kerr hole formation depends, of course, on the physics of angular momentum transport inside the progenitor 
star. The required spin could also be generated when the core mergers with a binary companion during common envelope phase. So 
perhaps one important distinction between a GRB and an ordinary supernovae is whether a black hole or a neutron star is formed in 
the aftermath. However, not all black hole formation events can lead to a GRB: if the minimum mass of a single star that leads to the 
formation of a black hole is as low as 25M Q , this would overproduce GRB by a large factor (see Izzard et al. 2004). 

5. THE AFTERGLOW APPEARANCE 

The interaction of a relativistic blast wave with the surrounding medium is described by the adiabatic Blandford & McKee (1976; 
hereafter BM) self-similar solution. The scaling laws that are appropriate for the burst interaction with a medium with particle 
density n oc r~ s have been described by Meszaros et al. (1998), Chevalier & Li 2000, Panaitescu & Kumar (2000), Wang et al. 
(2000), Ramirez-Ruiz et al. (2001), Dai & Lu (2002), and Dai & Wu (2003). 

For an adiabatic ultra-relativistic blast wave, the (isotropic equivalent) total energy is 

8irAT 2 r 3 - s c 2 

E= , (15) 

17 -4s ' 

where T is the bulk Lorentz factor of the shock front and r is the observed radius near the line of sight (BM). We assume the burst 
to be collimated with an initial half-angle 9 larger than 20°, and that lateral expansion is negligible during the relativistic phase. A 
distant observer receives a photon emitted along the line of sight towards the fireball center at a time (Chevalier & Li 2000) 

f = 4(4-Vv (16) 
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and so 

~(4-s)(n-4s)Et 



(17) 



2irAc 

Before the collision with the high-density shell, the shock front is expected to propagate through an n(r) = Ar~ 2 wind, where 

A = 3 x 10 35 M wr _5V~ 1 r 3 cm" 1 (18) 

for a spherically symmetric wind ejected at a constant speed. The relativistic expansion is gradually slowed down until the shock 
front encounters the wind density discontinuity at a radius r s h. A transition in the afterglow is then observed 

f sh w E52 >4, n A 35 days ( 1 9) 

after the burst. This of course only true for r s h < r^R (see equation 3). This encounter produces two new shock waves: a forward 
shock that moves into the wind thin shell discontinuity and a reverse shock that propagates into the relativistic ejecta. In the section 
that follows we examine the state of these shocks under various assumptions regarding the density discontinuity. 

5.1. Shock Conditions 

Consider a relativistic shell with a Lorentz factor r\ moving into the cold circumburst medium (CM). The interaction is described 
by two shocks: a forward shock propagating into the CM and a reverse shock propagating into the ejecta. There are four regions 
separated by the two shocks: the CM (1), the shocked CM (2), the shocked ejecta (3) and the unshocked ejecta (4). The CM is at rest 
relative to the observer (Fig. 5). Velocities and their corresponding Lorentz factors 7,, distances and time are measured relative to 
this frame. Thermodynamic quantities n,, p, and e, (particle number density, pressure, and internal energy density) are measured in 
the fluids' rest frames. The ISM and the unshocked shell are cold and, therefore, e\ = = 0. The shocked material is extremely hot 
and, therefore, px = ex/3 and pj = ej/3 (Sari & Piran 1995). 

For 77 = 74 » 1 the equations governing the shocks are (BM) 

— e * = 72 - 1 « 72; — = 4 72 + 3 «4 72 (20) 
ninipC 2 - n\ 

=73-1; - = 473 + 3. (21) 
nj,m p c z «4 

The approximations in equation (20) use only the fact that 74 3> 1 and therefore 72 S> 1 . No assumption is made about 73, the Lorentz 
factor of the motion of the shocked material in region (3) relative to the unshocked shell in region (4). 
Equality of pressures and velocities along the contact discontinuity yields 

«2=« 3 ; 7 3 R 4(- + -V (22) 

2 V72 74/ 

The solution for 72 depends only on two parameters, 77 and = rn/ni. The energy, pressure, and density also depend linearly on a 
third parameter, the external density n\ (Sari & Piran 1995). 

There are two simple limits of equations (20)-(22) in which the reverse shock is either Newtonian or ultra relativistic (the forward 
shock is always ultra relativistic if 77 3> 1 and ip > rj~ 2 ). If i] 2 » ip, the reverse shock is ultra relativistic (73 >• 1): 

73 = 72 = 73 = ~vr- (23) 

In this case almost all of the initial kinetic energy is converted by the shocks into internal energy (77 3> 73). The process therefore is 
over after a single passage of the reverse shock through the ejecta (Sari & Piran 1995). 
litp^>r] 2 the reverse shock is Newtonian (73 - 1 <C 1) and 

4t7 2 7/T 1 

73-1 « -y-=2e «1; 72 = 73 = ^7(1 -V^)- (24) 

The reverse shock converts only a fraction 7/ "C 1 of the kinetic energy into internal energy. It is too weak to slow down the ejecta 
effectively, and most of the initial energy is still kinetic energy when this shock reaches the inner edge of the ejecta. At this stage a 
rarefraction wave begins to propagate toward the contact discontinuity. This wave propagates at the speed of sound y/4pi /On^nip) 

and it reaches the contact discontinuity at t r = (3\/7 /4)Ar/\/ip/c, which is of the same order of magnitude as the shock crossing 
time. Here A is the width (in the observer's frame) of the relativistic shell. The rarefraction wave is then reflected from the contact 
discontinuity and a second, weaker, shock wave forms. A quasi-steady state slowing down solution forms after a few crossing like 
this (Sari & Piran 1995). Using momentum conservation, the total slowing down time can be estimated by ~ jAn4m p c/ p2- During 
this time the forward shock collects a fraction 7"' of the shell's rest mass. In contrary to the relativistic case, there are two relevant 
timescales now: the rarefraction or shock crossing time and the total slowing down timescale. 

In any realistic situation the CM is probably inhomogeneous, as in the stellar models described in the previous section. Consider 
a density jump by a factor a over a distance A. The forward shock propagates into the CM with a density n\ as before, and when 
it reaches the position where the CM density is an\ a new shock wave is reflected. This shock is reflected again off the shell. As 
discussed above, the reflections time is about t a ~ \/{4ca l l 2 ), and after these reflections, the Lorentz factor and hydrodynamical 
properties of the system are as if the CM were homogeneous with a density an\. If the reverse shock converts only a small fraction 
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of the kinetic energy into thermal energy, then the forward shock expands almost at the velocity of the previous blast wave. In the 
relativistic case, however, considerable deceleration may occur before t a unless ?a ~ r/y^pA/c > t a . 

The peak of synchrotron emission occurs at frequency v m k> q e j 2 B/(m e c) and the frequency integrated emissivity is given by 
e w <7Tcn e ^B 2 /Sir; where q e , m e , and j e are electron charge, mass, and thermal Lorentz factor, respectively, B 2 /8tt = eBii e j e m e is 
the magnetic field in the fluid rest frame (assumed to be amplified up to a fraction es by the processes in the shocked region and 
not determined directly by the field in the WR star), n e is the electron number density, and a T is the Thomson cross-section. The 
emission at low frequencies (y <C f m ) scales as v 1 ^ and at high frequencies (v 3> f m ) it scales as i/" ( '' _1) / 2 , where p w 2.5 is the 
power-law index for the energy distribution of electrons. The ratio of the peak synchrotron frequencies, before and after the forward 
shock has traversed a density jump of a factor a over a distance A, is 

5/2 



<?2,o 
<?2 



»2 
tl2,a 



-1/4 



(25) 



Therefore, we expect this emission to be a significant contribution to the long-wavelength flux 8 . The ratio of the observed flux after 
the forward shock has transversed through the density discontinuity, and the flux in the absence of the collision at a frequency much 
greater than the peak of the emission is given by 

where j a = 72,0/72 is a factor by which the Lorentz factor of the outer shell decreases as a result of the collision. The dependence 
of the observed flux on j a can be more rapid than the linear function considered above, depending on the temporal profile of the 
deceleration of the ejecta (see e.g. §5.2). The decrease in the observed emission from the forward shock due to the density jump at 
a frequency much greater than the peak of the emission is therefore given by a~* {p ~ l) w a~it, neglecting the enhanced energy losses 
from the shocks that may arise during the collision. This is clearly not the case at frequencies either below v m or for which there is 
significant emission from the reverse shock (typically a factor of 77 lower in frequency). At these frequencies, the observed flux is 
expected to increase. 



5.2. Afterglow Lightcurves 

We now generalize the above results to a spatially varying external density. The afterglow modeling used here is similar to that 
described in Salmonson (2003). The jet deceleration is calculated from the mass and energy conservation equations. The lateral 
spreading of the jet is neglected. The calculation of radiative losses includes synchrotron emission, and the synchrotron spectrum 
is taken to be a piecewise power law with the usual self-absorption, cooling, and injection break frequencies, calculated from the 
cooled electron distribution and magnetic field. The observed flux is obtained by integrating the jet emission over the equal arrival 
time surface. 

Strong temporal variations compared to the canonical power-law decay can be produced by changes in the circumburst density 
or by energy variations (e.g. Meszaros et al. 1998). For example, as shown in Ramirez-Ruiz, Merloni & Rees (2001), a shock 
wave that has been slowed by the surrounding medium could be caught up by subsequent shocks, increasing the shock wave energy. 
Alternatively, the shock wave's energy may have varied as the result of encountering an external medium of variable density. While 
here we concentrate on the case where the dominant variations are in density, the two alternatives can be distinguished by exploiting 
the fact that the flux at v > v c is insensitive to variations in the ambient density, whereas below the cooling frequency F v (v < v c ) oc n 1 / 2 
(Nakar et al. 2003; Heyl & Perna 2003). As the blast wave expands into the ambient matter, its kinetic energy is used to shock and 
heat the matter. Deceleration due to this starts in earnest when about half the initial energy is transferred to the shocked matter. The 
characteristic mass Md ec where this takes place is given by equation (2). This phase ends when so much mass shares the energy that 
the Lorentz factor drops to 1. This occurs at a mass scale M NR (see equation 3). 

In the case depicted in Fig. 3, for example, the stellar wind is not dense enough to slow down the ejecta to non-relativistic speeds 
before reaching r s h, so that we expect part of the blast wave evolution as we see it to take place outside the 1/r 2 density distribution. 
In a rarefied wind such as that illustrated in Fig. 4, the shock front will expand within the stellar wind until it reaches the density 
discontinuity at about 0.2 pc without significant deceleration. Over the typical time of observation of a GRB afterglow, the impact 
of the relativistic ejecta with the density enhancement will produce a clear feature in the observed emission. Fig. 6 and its caption 
summarize the predicted R-band afterglow lightcurve from hours until about a few years after the burst. One can see sharp break 
in the flux decay curve that coincides with a precipitous drop in Lorentz factor when the afterglow shock meets the shell. In this 
model, the isotropic energy of the ejecta is 5 x 10 53 erg. The characteristic synchrotron frequency is lower than that in the absence 
of the collision. This effect is responsible for the decrease in flux (at a fixed frequency) seen at about <~ 10 days. Further observable 
transitions are produced as the blast-wave plows deeper into the shocked-wind discontinuity. Jet effects are expected to become 
important at a time 

tj w 5(6» j /0.08) 4 (£/10 53 erg)(A/10 35 cm- 1 r 1 days, (27) 

just before the impact takes place. This makes the jet-break time not easily identifiable. 

The evolution of the apparent source size is shown in Fig. 7. The afterglow image is limited to a circle on the sky, whose size 
grows as f( 5 ~ s )/ 2 ( 4 ~ s ) (Granot & Loeb 2001). The assumption of a spherical flow may also serve as an adequate description of a jetted 

8 The above estimate assumes that the reverse shock is initially relativistic, which is likely to be the case when the density contrast of the density discontinuity is high 
(Dai & Lu 2002). 
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flow, at sufficiently early times before the jet break time, tj, when the Lorentz factor of the flow drops below the inverse of the jet 
opening angle (Rhoads 1997). In the case of a jet, the image at t > tj is expected to be different than in the spherical case (Oren, 
Nakar & Piran 2004; Granot, Ramirez-Ruiz & Loeb 2004) and will no longer be circular for observers who are situated away from 
the jet axis (Fig. 7). The image size depends decisively on the density profile of the ambient medium into which the GRB fireball 
propagates. The effects of the density discontinuity on the size of the afterglow image can be seen by comparing the solid and dashed 
curves in Fig. 7. This is not altogether surprising: the ejected material is decelerated by the external medium at a smaller radius than 
it would be in the absence of the density discontinuity. The calculations above demonstrate how the measurable properties of such 
afterglows depend on the nature of the progenitor star and the medium around it. 

6. CLUMP FORMATION AND AFTERGLOW VARIABILITY 

The smoothness or dumpiness of the swept-up shell may be responsible for small-scale variability observed in many gamma-ray 
burst afterglows. As illustrated in Fig. 4, the required degree of clumping in the Wolf-Rayet star wind itself does not seem plausible. 
On the other hand, as shown in Fig. 3, there is the possibility of clump development from the swept up RSG wind. As stated 
above, the shocked RSG wind is accelerated by the (less dense) shocked WR wind and is therefore subject to RT instabilities. The 
shell is initially thin and is therefore also subject to Vishniac instabilities. The two instabilities must then be coupled. The coupling 
mechanism appears to be a displacement of gas towards the interior of the bubble at the contact surface by the RT instability, which 
induces the formation of a valley in the outer shock, driving what is fundamentally a Vishniac instability. Garcfa-Segura et al. 
(1996b) found that clump formation is not efficient for the case of pure RT instabilities in thick shells and concluded that a necessary 
condition for clump formation to occur is that a shell be thin enough to allow Vishniac-like instabilities to drive transverse motions 
to form the clumps. Many RT fingers will then warp and break, forming knots. 

Once formed, the knots tend to dissipate. Several process tend to destroy the inhomogeneities embedded in a wind-blown bubble, 
namely thermal evaporation (Cowie & McKee 1977), ablation (Hartquist et al. 1986) and photo-evaporation (McKee 1986). If not 
inhibited by the presence of a magnetic field, heat conduction may deplete the clumps. The thermal evaporation time is 

n c r sw Ar 

f ev oc , (28) 

where T sw is the temperature of the hot shocked WR wind, Ar and « c are, respectively, the radii and densities of the clumps, and 
n sw is the shocked wind density. For n c > 10n sw , T sw ~4x 10 7 K and Ar < 0.01 pc one has f ev < 10 4 yr, so that, if conduction 
is unimpeded, the lifetime of the clumps may be shorter than the nebula lifetime. Erosion via ablation is far less effective, and the 
clumps survive if, as expected (see e.g. Hartquist et al. 1986), heat conduction is inefficient. The ablation time for a knot is 

n c Ar 

f a b« — • (29) 

^sw^rsg 

We estimate f ab > 10 2 f ev . The characteristic photo evaporation time of a clump at a distance r from a star emitting S ionizing photons 
per second is 

Ar 3 / 2 

f pe (r)cx« c -^ T y r r. (30) 

For the clump parameters considered above and S = 10 49 s" 1 one has f pe <~ 10 4 r pc yr. Small neutral condensations are thus likely to 
be photo-evaporated while larger ones survive. 

The jet will then encounter clumps - if they are not efficiently destroyed - with large variations in density. Presumably the clump 
is symmetric such that 

Ar,, =Ar ± = Ar clump . (31) 

Whether or not the ambient object collapses onto the shell is a key distinction that must be made in order to understand how the 
clump determines the time structure. We assume here the case of a "collapsible" object (such as a wind clump). The contribution to 
the peak duration from the time the jet takes to move through the clump Ar c i ump /(27 2 c), is negligible compared with the time the 
shell takes to engage the perpendicular size of the object. This engagement time is caused by the curvature of the jet. The curvature 
of the expanding jet prevents it from engaging the cloud instantaneously. Rather, the portion of the jet at 9 <~ 7" 1 requires a time 
r( 1 - cos 9)/ v longer to reach the cloud. At an angle 9 from the line of sight, the time to engage the object is 9 Ar c i ump / 2c. At a typical 
angle of 9 <~ 7" 1 , 

Af= 2Ar clump (32) 

C7 

The alternative is that the ambient source does not collapse, but produces photons on a scale of Ar c i ump (unlikely). The timescale 
in this case is determined from the light travel time across the overlap of the shell thickness and the ambient source thickness. If 
the clumps where distributed as a power-law of clump sizes, one might expect a power law distribution of observable peak widths. 
The cumulative effect of many small-scale density perturbations will tend to average out during the expansion history of the jet, 
and, therefore, the resultant afterglow profile should be relatively smooth. An upper limit on Ar c i ump is set by observed afterglow 
variability timescales (i.e. Af D b s ), where 

cAf obs < Ar dump < cA ^ bs7 , (33) 

Short time scale of afterglow oscillations provide interesting upper limits of Ar c i ump < 10(7/10) AU on the size of the clumps around 
GRB 011211 (Jakobsson et al. 2004) and GRB 030329 (Lipkin et al. 2004). These limits are lower than the fluctuation amplitudes 
seen on similar scales in the local interstellar medium (Wang & Loeb 1999), though they may reflect the length scale of comet-like 
clumps produced in ring nebulae surrounding massive stars. 
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7. DISCUSSION 

It is evident from the above discussion that the environment of a massive star at the time of its death is a very rich one. Even in the 
simplest case of a wind whose properties do not vary over the life of the star, complex behaviour with multiple possible transitions 
in the observable part of the afterglow lifetime may be seen. The eventual resulting afterglow lightcurve depends fairly strongly on 
the properties of the system, especially the mass-loss rate of the star and the ambient density. This has a good and a bad side. On 
the negative side, it implies that one can not be too specific about the times at which we expect to see transitions in the observed 
emission. On the positive side, if and when we do see these transitions, they can be fairly constraining on the properties of the stellar 
progenitor. 

Considering the radial range relevant to GRBs, the absence of the expected 1/r 2 density structure in many bursts is not surprising. 
On the other hand, the reason for low densities remains unclear. A wind termination shock may resolve this matter although very 
special conditions are necessary to bring this about. Because the density in a shocked wind is higher than that in a free wind at the 
same radial distance, the low-density requirement is not alleviated by appealing to a shocked WR wind. A low stellar mass wind loss, 
a faster wind velocity or a low metallicity may help, although it is not clear whether realistic assumptions can provide the required 
low densities (Wijers 2001; Chevalier et al. 2004). The following point should be emphasized here. The fireball model used to infer 
the circumburst density is highly simplified. For example, the wind density is assumed to follow a pure 1/r 2 law and to be free of 
inhomogeneities, the expanding jet surface is assumed uniform with no internal (density, velocity) gradients and the fraction of the 
explosion energy in the post shock magnetic field is assumed to be constant (Panaitescu & Kumar 2002). We do not know whether 
a description of the afterglow data is possible with a wide variety of underlying assumptions and whether that would substantially 
change the parameters inferred. Still, the low densities inferred by afterglow observations are thus problematic for the collapsar 
model. 

The task of finding useful progenitor diagnostics is simplified if the pre-burst evolution leads to a significantly enhanced gas 
density in the immediate neighborhood of the burst. The detection of spectral signatures associated with the GRB environment 
would provide important clues about the triggering mechanism and the progenitor. A special case is that of GRB 021004, where 
lines of highly ionized species, blueshifted relative to the host galaxy, have been attributed to a Wolf-Rayet stellar wind (Mirabal et al. 
2003; Schaefer et al. 2003). Stars interact with the surrounding interstellar medium, both through their ionizing radiation and through 
mass, momentum and energy transfer in their winds. Mass loss leads to recycling of matter into the interstellar medium, often with 
chemical enrichment. Mass loss is a significant effect in the evolution of massive stars; in particular, for stars that enter WR stages 
(e.g. Ramirez-Ruiz et al. 2001). WR stars follow all or part of the sequence WNL, WNE, WC and WO, which corresponds to a 
progression in the exposure of nuclear products (CNO equilibrium with H present, CNO equilibrium without H, early visibility of 
the products of the 3a reaction and then a growing (C+0)/He ratio respectively). Mass loss drastically influences stellar yields. At 
low Z there is a high production of a-nuclei, while at higher Z most of the He and C produced is ejected in stellar winds and escapes 
further nuclear processing. This characteristic offers a direct observational test of which stars are likely to produce a GRB. 

Finally, The total energy observed in 7-rays from GRBs whose redshift has been determined is diverse (e.g. Soderberg et al. 2004). 
One appealing aspect of a massive star progenitor is that the great variety of stellar parameters can probably explain this diversity. 
Given the need for a large helium core mass in progenitors, the burst formation may be favored not only by rapid rotation but also 
by low metallicity. Larger mass helium cores might have more energetic jets, but it is unclear whether they can be expected to be 
accelerated to large Lorentz factors (MacFadyen, Woosley & Heger 2001; Ramirez-Ruiz, Celotti & Rees 2002). Many massive stars 
may produce supernovae by forming neutron stars in spherically symmetric explosions, but some may fail neutrino energy deposition, 
forming a black hole in the center of the star and possibly a GRB. One expects various outcomes ranging from GRBs with large 
energies and durations, to asymmetric, energetic supernovae with weak GRBs (Totani 2003; Granot & Ramirez-Ruiz 2004). The 
medium surrounding a GRB would provide a natural test to distinguish between different stellar explosions. 

We thank M. J. Rees, J. Granot, R. Chevalier, S. Woosley and N. Langer for helpful discussions. We are particularly grateful to 
the anonymous referee for helpful comments. This work is supported by NASA through a Chandra Postdoctoral Fellowship award 
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FIG. 2. — The formation of shocks by a stellar driven wind interacting with the surrounding interstellar medium. The interaction leads to a driven wave composed 
of shock compressed wind, a contact surface, and the swept-up ISM. Fast stellar wind matter (zone [a]) enters the shock where it is compressed by a factor of 
~ 4, and heated; the accumulated shock wind material is in zone (b). The entire region between the two shocks is nearly isobaric, so if material cools, it becomes 
compressed and resides in the contact surface, which forms a boundary between the shocked wind and the shocked, swept-up ISM (zone [c]). The undisturbed ISM 
is in zone (d). 
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FIG. 4. — The state of the nearby circumstellar medium around a 60 Mq rotating massive star. During the last stages within its WR lifetime, the wind velocity 
increases. The fast wind collides with early ejecta before reaching the LBV shell. The reader is refer to Fig. 2 of Garci'a-Segura et al. (1996b) for the stellar 
mass-loss rates and wind velocities as a function of time. The grid has 800 X 1 80 zones, with a radial extent of 0.5pc and an angular extent of 90° . The inner-most 
radial zone lies at 0.0125 pc. Left panel: Logarithm of the circumstellar density in units of g era" 3 . Right Panel: Density and cumulative mass as a function of radius 
along the polar axis. 
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FIG. 5 . — Diagram illustrating a relativistic shock system. Basic system consisting of a shocked fluid encountering matter at rest. Quantities for the system are the 
unshocked ejecta, the reverse and forward shocks, and the external matter at rest. The dashed line schematically shows the properties of the relativistic ejecta in the 
absence of an external medium. 
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FIG. 6. — The effect of the impact of a relativistic jet with the wind density discontinuity on the R-band afterglow. The shock front expands within the .v = 2 
stellar wind until it reaches the high-density shell at a distance r s h ss 0.2 pc (see Fig. 3). The shock transit of the massive shell causes a rapid decline in Lorentz 
factor and a corresponding decline in flux (at a fixed frequency). The remaining evolution of the shock is non-relativistic. At the time of the collision the relativistic 
shell Lorentz factor is 7 ~ 4 for E = 5 X 10 53 ergs. In this simulation the jet opening angle is 9j = 5° and it is viewed at o t> s = 3° from the jet axis. The afterglow 
emission is calculated in the adiabatic regime. The collision model takes into account the fireball geometrical curvature when calculating the photon arrival time and 
relativistic boosting. 
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FIG. 7. — Evolution of the afterglow source size as a function of time for a sharp-edged, homogeneous jet seen at 9^ = 0.69j = 3°. This one dimensional, 
hydrodynamic model takes on the density profile of the stellar wind in the polar direction. The jet deceleration is calculated from the mass and energy conservation 
equations. The lateral spreading of the jet is neglected. The effect of the impact of a relativistic jet with the wind density discontinuity (Fig. 3) can be seen by 
comparing the solid and dashed curves. For simplicity, the density profile seen in Fig. 3 has been extrapolated to large radii. Substantial deviations in the density 
profile are, however, expected to occur with the inclusion of the wind termination shock region (see e.g. Fig. 2 in Garcla-Segura et al. 1996a). The inset panel shows 
the afterglow source diameter. The x- and y-axes are shown in their true scale and measured in light-days. 



